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ABSTRACT 

MSH 15-56 (G326.3-1.8) is a composite supernova remnant (SNR) that consists of an SNR shell and 
a displaced pulsar wind nebula (PWN) in the radio. We present XMM-Newton and Chandra X-ray 
observations of the remnant that reveal a compact source at the tip of the radio PWN and complex 
structures that provide evidence for mixing of the supernova (SN) ejecta with PWN material following 
a reverse shock interaction. The X-ray spectra are well fitted by a non-thermal power-law model whose 
photon index steepens with distance from the presumed pulsar, and a thermal component with an 
average temperature of 0.55 keV. The enhanced abundances of silicon and sulfur in some regions, and 
the similar temperature and ionization timescale, suggest that much of the X-ray emission can be 
attributed to SN ejecta that have either been heated by the reverse shock or swept up by the PWN. 

We find one region with a lower temperature of 0.3 keV that appears to be in ionization equilibrium. 
Assuming the Sedov model, we derive a number of SNR properties, including an age of 16,500 yr. 
Modeling of the 7 -ray emission detected by Fermi shows that the emission may originate from the 
reverse shock-crushed PWN. 


1. INTRODUCTION 

Composite supernova remnants (SNRs) are those con- 
sisting of a central pulsar that produces a wind of 
synchrotron-emitting relativistic particles, and a super- 
nova (SN) blast wave that expands into the surrounding 
interstellar medium (ISM). The interactions between the 
SNR and the evolving pulsar wind nebula (PWN) can re- 
veal a wealth of information about the SNR properties, 
the SN ejecta, the central pulsar, and the spectrum of 
particles injected into the PWN. At the late stages of a 
composite SNR’s evolution, the SN reverse shock crushes 
the PWN, resulting in complex filamentary structures 
and mixing of the PWN material with ejecta gas (e.g. 
Blondin et al. 2001; van der Swaluw et al. 2004). In 
cases where the reverse shock interacts with the PWN 
asymmetrically, either due to the pulsar’s motion or a 
density gradient in the ambient ISM, the PWN can be 
swept away from the pulsar, resulting in a relic PWN 
(e.g. van der Swaluw et al. 2004). Due to the burn-off 
of high-energy particles, and the fact that new particles 
are no longer being injected into the relic PWN, their 
emission is usually observed at radio wavelengths. As in 
the case of G327. 1-1.1 (Temim et al. 2009), the pulsar 
and its newly forming PWN are displaced from the relic 
nebula. The interaction between the reverse shock and 
the PWN may also play an important role in the origin 
of gamma-ray emission in composite SNRs (e.g. Slane 
et al. 2012). The increase in the magnetic field due to 
such an interaction produces an excess of low energy par- 
ticles that may give rise to gamma-ray emission through 
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inverse Compton scattering. 

With its first detection at radio wavelengths, the SNR 
MSH 15-56 (G326.3-1.8) became the prototype for the 
composite class of SNRs (Mills et al. 1961; Milne et al. 
1979). A lower limit on its distance of 3.1 kpc has been 
established by the H I absorption profile of Goss et al. 
(1972), while Rosado et al. (1996) found a distance of 
4.1 kpc through Ha velocity measurements, the distance 
we adopt in this paper. Higher-resolution radio obser- 
vations taken with the Molonglo Observatory Synthesis 
Telescope (MOST) (Whiteoak & Green 1996) and the 
Australia Telescope Compact Array (ATCA) showed a 
symmetric SNR shell, and plerionic component that is 
displaced from the geometric center of the SNR (Green 
et al. 1997; Dickel et al. 2000), reminiscent of a relic 
nebula that has been disrupted by the reverse shock. 
The radio emission is non-thermal, with a spectral index 
of a = 0.34 for the non-thermal shell emission, where 
S v oc v~ a , and 0.18 for the PWN component (Dickel 
et al. 2000). The PWN component is highly polarized, 
with a luminosity of L 10 7_ 10 iihz ~ 5 x 10 34 ergs _1 , mak- 
ing the radio PWN in MSH 15-56 the third most lumi- 
nous after the Crab and G328.4+0.2 (Dickel et al. 2000). 
A search for the pulsations from the pulsar that powers 
the nebula have been unsuccessful (Kaspi et al. 1996). 

Optical Ha filaments were observed on the face of the 
SNR by van den Bergh (1979), and appear to correlate 
spatially with the SNR shell. MSH 15-56 was also de- 
tected at X-ray wavelengths by ROSAT (Kassim et al. 
1993) and ASCA (Plucinsky 1998). The X-ray morphol- 
ogy shows a complete shell that spatially correlates with 
the radio SNR. The emission appears to be clumpy, with 
an enhancement in the regions near the PWN. Spectral 
fits to the ROSAT data showed that the global X-ray 
emission from 0. 1-2.4 keV can be fit with a column den- 
sity of Nh = 8.9 x 10 21 cm -2 and a thermal model with 
a temperature of 0.56 keV (Kassim et al. 1993). 

In this paper, we report the analysis of the XMM- 
Newton and Chandra observations of MSH 15-56 that 
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provide insight into the nature of the X-ray thermal emis- 
sion and the evolutionary state of the remnant. We 
also model the Fermi Gamma-ray Space Telescope Large 
Area Telescope {Fermi- LAT) gamma-ray emission that 
appears to be associated with the SNR and discuss its 
possible origins. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. X-ray 
2.1.1. XMM-Newton 

XMM-Newton observations of MSH 15-56 were car- 
ried out on 2004 August 11 with the MOS and pn cam- 
eras, under the observation ID 020427010. The obser- 
vations were taken in two different exposures for a total 
exposure time of 56 ks. However, the observations were 
significantly affected by flares, so after the standard pro- 
cessing with the XMM-SAS software, version 11.0.0, the 
resulting exposure times for the MOS 1, MOS2, and the 
pn instruments were 29.9 ks, 29.7 ks, and 23.0 ks, respec- 
tively. 

The cleaned data from all three instruments were com- 
bined to create exposure-corrected images in three differ- 
ent energy bands using the XMM-SAS images script. 7 
The image is shown in Figure 1. The spatial binning size 
for the images was set to 5". For the purpose of spectral 
extraction, data from both exposures of the two MOS in- 
struments were merged into a single MOS dataset using 
the merging recipe from the SAS website (website here) . 
In the same way, the two pn exposures were merged into 
a single pn dataset. The spectra were then extracted 
from the merged MOS and pn files using SAS version 
11.0.0, using the source and background regions shown 
in Figure 2a, and fitted simultaneously. The spectra were 
grouped to include at least 30 counts in each bin, back- 
ground subtracted, and fitted using the CIAO 4.3 Sherpa 
software. All spectra were fitted with a two-component 
model consisting of a power-law and the XSVNEI ther- 
mal component. The XSPEC model tbabs was used for 
the absorption along the line of sight. We found that 
the two-component model significantly improved the y 2 
value for all regions. 


2.1.2. Chandra 

MSH 15-56 was observed with the Chandra X-ray 
Observatory’s Advanced CCD imaging Spectrometer, 
ACIS-S. Observations were carried out on 2001 August 
18 for a total exposure time of 56.8 ks, under the observa- 
tion ID 1965. The standard cleaning and data reduction 
were performed using CIAO version 4.3. Since the emis- 
sion from the SNR fills the entire held of view of the ACIS 
detector, we extracted background spectra from the 
blank-sky hies produced by the ACIS calibration team 
(http: / / cxc.harvard.edu/ciao/threads/acisbackground/). 
A background spectrum extracted from the blank-sky 
data using the same region was subtracted from each 
spectrum, and the residual data were binned and 
htted using the CIAO 4.3 Sherpa software. The same 
two-component spectral model as for the XMM data 
was used to ht the Chandra spectra. 

7 http:/ /xmm.esac.esa.int/sas/current/documentation/threads/ 
epic_merging.shtml 



Fig. 1. — (a) MOST 843 MHz radio image of MSH 15-56 
(Whiteoak & Green 1996) in green, and the 0.1-2. 4 keV ROSAT X- 
ray image (Kassim et al. 1993) in violet. The region in the white 
dashed rectangle is shown in panel (b). (b) A combined three 

color XMM-Newton X-ray image (MOSl/MOS2/pn) with 0. 5-1.0 
keV emission shown in red, 1. 0-3.0 keV image in green, and 3. 0-9.0 
keV image in blue. MOST radio contours are overlaid in white. 
The region in the white dashed rectangle is shown in panel (c). 
(c) Chandra X-ray 0.3-10.0 keV image of the PWN in MSH 15-56 
resolving the trail behind a point source and an arc of emission 
tracing the edge of the radio PWN. The MOST radio contours are 
overlaid in white. 
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Fig. 2. — (a) A combined three color XMM-Newton X-ray image (MOSl/MOS2/pn) with 0. 5-1.0 keV emission shown in red, 1. 0-3.0 keV 
image in green, and 3. 0-9.0 keV image in blue. The spectral extraction regions are overlaid in white, along with the dashed background 
region, (b) Chandra X-ray 0.3-10.0 image with the spectral extraction regions overlaid in white. The figures are on the same spatial scale 
as those in Figure 1. 


2.2. 7 -Ray 

In this work, 51 months of data from Fermi - LAT (from 
August 2008 until November 2012) were analyzed. Only 
events belonging to the Pass 7 V6 Source class, which re- 
duces the residual background rate as explained in detail 
in Ackermann et al. (2012), were selected for this study. 
The updated instrument response functions (IRFs) called 
“Pass7 version 6” , developed using in-flight data (Rando 
et al. 2009; Ackermann et al. 2012), were used. Addition- 
ally, only events coming from zenith angles smaller than 
100° were selected so as to reduce the contribution from 
terrestrial albedo 7-rays (Abdo et al. 2009b). This study 
included data from a circular region in the sky centered 
on the position of MSH 15-56, with radius 20°, analyzed 
using the Fermi Science Tools v9r23pl 8 . 

The maximum likelihood fitting technique was em- 
ployed to study the morphological and spectral char- 
acteristics of the region (Mattox et al. 1996). The 
emission models used in gtlike included a Galactic dif- 
fuse component resulting from the interactions of cos- 
mic rays with the ISM and photons (using the mapeube 
file gal_2yearp7v6_v0_trim. f its), and an isotropic one 
that accounts for the extragalactic diffuse and residual 
backgrounds (modeled using the iso_p7v6source.txt 
table) . 

In order to understand the spatial characteristics of 
the 7-ray emission in the field of MSH 15-56, data in the 
energy range 2 to 200 GeV, and converted in the front 
section, were used. The 68% containment radius angle 
for normal incidence /roni-selected photons in this en- 
ergy band is < 0.4°. Galactic and isotropic backgrounds 
were modeled and test statistic (TS) maps were con- 
structed using gttsmap to allow for detection significance 
estimates, and to best evaluate the position and possible 
extent of the source. The test statistic is the logarithmic 


ratio of the likelihood of a point source being at a given 
position in a grid, to the likelihood of the model without 
the additional source, 2 log(L ps / L nu u). 

The study of the spectral energy distribution (SED) 
characteristics of the source associated with MSH 15-56, 
was performed combining events converted in both front 
and back sections, and in the energy range 0.2-204.8 GeV. 
The lower energy bound was selected to avoid the rapidly 
changing effective area of the instrument at low energies, 
and because of the large uncertainty below 0.2 GeV re- 
lated to the Galactic diffuse model used, gtlike is used to 
model the flux at each energy bin and estimate, through 
the maximum likelihood technique, the best-fit parame- 
ters. Background sources from the 24-month Fermi LAT 
Second Source Catalog (Nolan et al. 2012) 9 were included 
in the model likelihood fits. The systematic uncertainties 
of the effective area, for the IRF used, are energy depen- 
dent: 10% at 100 MeV, decreasing to 5% at 560 MeV, 
and increasing to 20% at 10 GeV (Abdo et al. 2009a; 
Nolan et al. 2012, and references therein). As an addi- 
tion to the statistical uncertainties associated with the 
likelihood approach, and the systematic errors related to 
the IRFs, the uncertainty of the underlying Galactic dif- 
fuse level was considered. This source of uncertainty was 
included in the evaluation of the systematic uncertainties 
by artificially varying the normalization of the Galactic 
background by ±6% from the best-fit value at each en- 
ergy bin, similarly to the analysis used in Castro & Slane 
( 2010 ). 

3. ANALYSIS AND RESULTS 
3.1. Multi- Wavelength Morphology 

SNR MSH 15-56 is the prototypical example of a com- 
posite SNR. Figure la shows the MOST 843 MHz radio 
image of MSH 15-56 (Whiteoak & Green 1996) in green, 


8 The Science Tools package and related documentation 

are distributed by the Fermi Science Support Center at 
http:/ /fermi. gsfc.nasa.gov/ssc 


9 The data for the 1873 sources in the Fermi LAT Second Source 
Catalog was made available by the Fermi Science Support Center 
at http: / / fermi.gsfc.nasa.gov/ssc/data/ access /lat / 2yr .catalog/ 
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and the 0. 1-2.4 keV ROSAT X-ray image (Kassim et al. 
1993) in violet. The radio morphology consists of a sym- 
metric filamentary shell, with a bright PWN component 
offset to the SW from the geometric center. The shell is 
roughly circular, and approximately 18' in radius, which 
corresponds to ~ 21.5 pc for a distance of 4.1 kpc. The 
shell emission is brightest in the NW, where the X-ray 
emission appears to be fainter. The shape of the PWN 
is irregular and suggests that the nebula has likely inter- 
acted with the SN reverse shock. It is elongated along 
the SE/NW direction, with approximate dimensions of 
516 along the long and 314 along the short axis. Dickel 
et al. (2000) analyzed the high resolution ATCA images 
of MSH 15-5 6 that revealed more detailed structures of 
the PWN, consisting of parallel ridges aligned with the 
direction of the long axis. The integrated flux densities 
at 1 GhZ are 114 Jy and 26 Jy for the shell and PWN 
components, respectively. 

The X-ray emission is patchy is appearance and fills the 
entire SNR shell. The ROSAT soft X-ray emission (0.1- 
2.4 keV) is seen in violet in Figure la. There are localized 
regions of enhanced emission, but they do not correlate 
with any radio structures. The emission in the whole 
southeastern half of the remnant and around the PWN 
is enhanced, but there is no spatial correlation between 
the diffuse X-ray emission and the radio PWN. Kassim 
et al. (1993) found that the X-ray emission is mostly 
thermal, with a temperature of 0.56 keV. 

The XMM-Newton observations are shown in Figure 
lb, with the 0. 5-1.0 keV emission shown in red, 1.0- 
3.0 keV image in green, and 3. 0-9.0 keV image in blue. 
The MOST radio contours are overlaid in white. The 
XMM field of view includes the region around the ra- 
dio PWN, and does not cover the entire shell. There 
is a point source located at the very tip of the radio 
PWN, and a trail of hard X-ray emission immediately 
behind the point source, in the direction of the radio 
PWN. The western emission appears softer, while the 
emission around the northern part of the PWN appears 
more green in the three-color image. There is a patch of 
very soft emission SE of the PWN, near the edge of the 
radio shell contour. 

The dashed white rectangle in Figure lb represents 
the Chandra field of view shown in Figure lc. The high- 
resolution Chandra image shows the detailed structure of 
the X-ray PWN. A point source located at the tip of the 
radio PWN, at the position 5 h 52 m 12.5^, -56°18'58"0. The 
trail of X-ray emission behind the point source has a bow- 
shock-like morphology, suggestive of a PWN produced 
by a high-velocity pulsar. Such cometary morphologies 
are observed in the Mouse, Geminga, G327. 1-1.1, and a 
number of other PWNe (e.g. Kargaltsev & Pavlov 2008). 
Similarly, there is an arc of X-ray emission tracing the 
southern edge of the contour of the radio PWN. Diffuse 
X-ray emission is present west of the point source, but 
this emission does not appear to correlate with any radio 
structures. 

3.2. X-Ray Spectroscopy 

X-ray spectra were extracted from the XMM and 
Chandra data from the regions shown in Figure 2. For 
the XMM data, we selected eight different regions that 
correspond either to specific structures identified in the 
Chandra image, or to regions that show differences in col- 


ors in the three-color XMM image. The selected spec- 
tral regions in Figure 2a are the following; (a) the arc 
structure that traces the edge of the radio PWN, (b) the 
cometary X-ray trail behind the point source, (c) the cen- 
tral region of the diffuse X-ray emission, (d) the softest 
X-ray patch that appears red in the XMM three-color 
image, (e) the northern region that appears most green 
in the XMM three-color image, (f, g) the diffuse X-ray 
emission east of the PWN, and (h) the larger, more dif- 
fuse arc of X-ray emission visible in the Chandra image. 
A region outside of the SNR’s radio contours was se- 
lected for the background exctraction, and this region is 
indicated by the dashed white ellipse in Figure 2a. The 
spectral extraction regions used for the Chandra data are 
shown in Figure 2b. The regions correspond to distinct 
structures; (a) the cometary X-ray trail, (b) the X-ray 
point-like source at the tip of the radio PWN, (c) the 
emission east of the PWN, and (d) the arc tracing the 
edge of the radio PWN. The background for the Chandra 
data was extracted from the blank-sky data, as described 
in Section 2. The extracted XMM and Chandra spectra 
are shown in Figures 3 and 4, respectively. 

We first attempted to fit the X-ray spectra with a sin- 
gle XSVNEI model, but found that a high-energy excess 
required an additional power-law component for most re- 
gions. We therefore fitted all the spectra with a two- 
component model, a power-law plus a thermal XSVNEI 
model, using the XSPEC absorption model TBABS. The 
point source spectrum was the only one that was fit- 
ted with a power-law component only. The absorption 
along the line of sight was linked and simultaneously 
fitted for all the XMM regions. The best fit value of 
Nh = 0.51 x 10 22 cm -3 was then used as a frozen param- 
eter for the fitting of the Chandra spectra. The abun- 
dances of silicon and sulfur were clearly enhanced in re- 
gion (e) of Figure 2a, the northernmost region in the 
image, so we thawed the silicon and sulfur abundances 
for all the fits in order to determine if the there is a possi- 
ble enhancement. The best fit parameters for the XMM 
and Chandra fits are listed in Tables 1 and 2. 

Based on the best fit parameters, we make the follow- 
ing observations: 

(1) All regions besides the soft X-ray patch (region 
(d) of Figure 2a) are fitted with approximately the same 
thermal temperature of around 0.5-0. 6 keV and the same 
ionization timescale of approximately 2 x 10 11 scm -3 . 

(2) Region (e) of Figure 2a shows enhanced abundances 
of silicon and sulfur that are up to 3 or 4 times the solar 
value, providing clear evidence that the emission in this 
region comes from SN ejecta. The arc feature in Figure 
2b also shows enhanced abundances, consistent with the 
interpretation that this emission comes from SN ejecta 
that have been swept up by the PWN. The abundances 
in other regions are not well constrained, so we cannot 
discount the possibility that the emission from other re- 
gions also comes from SN ejecta. 

(3) Region (d) is the only region for which the best-fit 
parameters of the thermal component significantly differ 
from the others. The best-fit temperature is 0.3 keV 
and the ionization timescale is significantly longer, on 
the order of 10 13 s cm -3 . 

(4) The photon index steepens with distance from the 
point source, varying from F = 2.3 in the cometary PWN 
structure to T = 4.5 in the outer regions. Similar steep- 
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Fig. 3. — XMM-N ewton X-ray spectra extracted from regions (a)-(h) in Figure 2a. The best-fit models are shown as solid red lines. The 
individual components are the XSVNEI thermal model (dotted red line) and a power- law model (dashed red line). 


TABLE 1 

XMM SPECTRAL FITTING RESULTS 


REGION: 

a 

b 

c 

d 

N H (10 22 cm- 2 ) 

n 51 +0-002 
u - D1 -o.ooi 




Photon Index (r) 

o qfi+ 0.09 

z.ou_025 

2 rn+°- 12 

q 1 o+O.OT 
O. 1 O- 0.07 

4 47 + 0-32 
4 - 4 '- 0.30 

Amplitude 

(l-ltoi) X 10- 4 

(7.5tS;g) X 10- 5 

(8.2++ x 10- 4 

(i-4±B:i) x io- 4 

Fi (erg/cm 2 /s) 

5.4 x 10 -13 

4.2 x 10 -13 

4.5 x IO' 12 

1.7 x IO" 12 

fcT(keV) 

0 5fi+°- 20 

U.O D _oo2 

0 62”*" 0,30 
U.OZ-o 20 

0.58l° 0 '“i 

0 OQ+ 0,04 
u - zy -0.0l 

Silicon 

1.6±?‘? 

A e+ 3.8 
4 *°- 2.3 

i-stSl 

1 7+O.8 
4 ‘ * —0.8 

Sulfur 

6-6±|;i 

0.0+ 3 

1 q+0.7 
0.8 

l.ot 4 ; 9 

r(s cm -3 ) 

(2.3+?+ x 10 11 

(i-it?:?) x io 11 

(1.8±j+) x 10 11 

(3.5±3.- 4 ) x 10 13 

Normalization 

(i-s±8:l) x io- 4 

(3.9+3) x 10~ 5 

(2.0±°J) x 10“ 3 

(2.2±g;|) x IO” 3 

F2 (erg/cm 2 /s) 

4.7 x IO -13 

1.8 x 10“ 13 

7.2 x 10~ 12 

4.3 x IO -12 

Reduced \ 2 Statistic 

0.7 

0.7 

1.1 

1.0 

REGION: 

e 

f 

g 

h 

N H (10 22 cm — 2 ) 
Photon Index (r) 

3-80 ±g;}^ 

9 qn+ 0.05 

Z.OU- 0.26 

4 54 + 1 - 98 
^•°^-0.20 

bO 

bo 

to • 
1 + • 
° 9 ’ 
o o 
-4 00 

Amplitude 

(8.4+3) x 10“ 4 

(i-7±g;|) x io- 4 

(2.5± < j ) §) x 10" 4 

(4.4t°; 4 ) x 10- 4 

Fl (erg/cm 2 /s) 

6.4 x IO” 12 

8.6 x 10“ 13 

3.3 x 10~ 12 

7.0 x 10^ 13 

fcT(keV) 

0 66 +0 ' 01 
u.oo_o 07 

0 51+ 0 - 10 
U. 01-0 01 

n 5 i +°- 14 
u - 0i -0.01 

0 54 + 0 - 01 
u,o ^-0.01 

Silicon 

2 9 +0 ' 7 
Z-y -0.2 

i 9+O.2 
1,z — 0.6 

0 Q+°- 2 
U,y -0.5 

1 3 +0 - 4 

Sulfur 

o rr+ 0.2 

1.0+ 1 ' 4 

0.8 +1 '° 

0.0+ 2 ' 1 

r(s cm -3 ) 

(1.6toi) x 10 11 

(1.9±g'|) x 10 11 

(2.1+ 2 ) x 10 11 

(2.0+ 4 ) x 10 11 

N ormalizat ion 

(2.2+J®) x IO" 3 

(7.0±g*J) x 10~ 4 

(1-2++ x IO" 3 

(1.4++ x IO" 3 

F2 (erg/cm 2 /s) 

8.7 x 10~ 12 

2.5 x 10“ 12 

4.1 x 10~ 12 

1.7 x IO’ 12 

Reduced \ 2 Statistic 

0.9 

0.8 

1.1 

0.8 


Note. — The listed uncertainties are 1-sigma statistical uncertainties. The apertures used in the 
spectral extractions are shown in Figure 2a, and the spectra and fits are shown in Figure 3. The 
fluxes were calculated in the 0.3-10 keV band. The normalization of the thermal models is equal to 
10 — 14 n e nif V/4:7rd 2 cm -5 , where V is the volume of the emitting region and d is the distance to the 
SNR. 


6 


TEMIM ET AL. 


0.1 

> 0.01 

% 0.001 
o 

0.0001 

3 
0 
-3 
0.1 

0.01 

0.001 

0.0001 

3 

(0 

o) 0 

in 

-3 


03 

E 

O) 

in 


> 

0 

-X. 


o 



1 


10 


10 


Energy (keV) 


Energy (keV) 


Fig. 4. — Chandra X-ray spectra extracted from regions in Figure 2b. The best-fit models are shown as solid red lines. The individual 
components are the XSVNEI thermal model (dotted red line) and a power-law model (dashed red line). 


TABLE 2 

CHANDRA SPECTRAL FITTING RESULTS 


REGION: 

a 

b 

c 

d 

N H (10 22 cm- 2 ) 

0.51 




Photon Index (T) 

1 84”^ 0 ' 13 

I 41+0.24 

2 n^+O-H 
z ’ uo -0.08 

1 60 +0 ' 12 
± - ou -0.13 

Amplitude 

(5.2 ±°'g) x 10- 5 

(5.6^j' 3 ) x 10-® 

(2.1±g.2) x 10 -4 

(6.6±g !) x 10- 5 

Fi(erg/cm 2 /s) 

3.2 x 10 -13 

5.2 x 10 -14 

1.1 x 10 -12 

5.0 x 10- 13 

fcT(keV) 

0 49 +0 10 

u.^y_ 0 .i3 


O.DO_ 0> 02 

0 62”*" 0 ' 04 
u - dz -0.04 

Silicon 

o q+ 4.6 
c5 - y -2.6 


1 q+0.2 

i ' U -0.3 

o o+0.6 
■^•0-0.6 

Sulfur 

o.ot 13 


l.ot 1 - 4 

5-6±l;g 

r(s cm -3 ) 

(2.5±i. 3 ) x 10 11 


(l.8t°; 3 ) x 10 11 

(l.7l° 0 ; 4 4 ) x 10 11 

N ormalizat ion 

(4.0±5 |) X 10~ 5 


(5.9±g;l) x 10~ 4 

(l-4±ol) x 10~ 4 

F 2 (erg/cm 2 /s) 

1.3 x 10 -13 


2.1 x 10~ 12 

5.4 x 10~ 13 

Reduced \ 2 Statistic 

0.7 

0.9 

0.9 

0.6 


Note. — The listed uncertainties are 1-sigma statistical uncertainties. The apertures used in the 
spectral extractions are shown in Figure 2b, and the spectra and fits are shown in Figure 4. The 
fluxes were calculated in the 0.3-10 keV band. The normalization of the thermal models is equal to 
10 — 14: n e riH V/And? cm -5 , where V is the volume of the emitting region and d is the distance to the 
SNR. 
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ening with distance from the pulsar has been observed 
in other PWNe, including 3C 58, G21. 5-0.9, G327.1-1.1, 
and G54. 1+0.3 (Slane et al. 2000, 2004; Temim et al. 
2009, 2010). The steepening can be explained qualita- 
tively by synchrotron aging of particles as they travel 
from the pulsar towards the outer nebula (Kennel & 
Coroniti 1984). The large photon index of T > 3.5 that 
is seen in the outermost regions labeled (e), (d), and (g) 
in Figure 2a is somewhat steeper than expected. We 
suggest that an underlying soft X-ray component from 
swept-up ISM is contributing to the X-ray spectrum and 
steepening the apparent photon indices of the power-law 
component (see Section 6). Even though this additional 
soft X-ray component cannot be constrained by the cur- 
rent data, our simulations show that the presence of such 
a component with a temperature of ~ 0.3 keV and a sur- 
face brightness three times lower than in region (d) would 
raise the temperature of the presumed ejecta component, 
and flatten the global photon index of the power-law 
component from T = 3.0 ± 0.1 to T = 2.4 ± 0.1. 

The total luminosity of the X-ray emission is Lx (0.3 — 
10 keV) = 4.0 x 10 34 erg s _1 . Using a value of E = 
4 x 10 36 erg s _1 (see Section 7), this results in an X-ray 
radiation efficiency of 0.01. While this value is somewhat 
higher than average, it falls well within the scatter in 
the Lx/E relationship observed for other pulsar (e.g. 
Possenti et al. 2002). 

(5) The spectrum of the point source is fitted with a 
power-law with a photon index of T = 1.4 ± 0.2, con- 
sistent with the value expected for a pulsar. The unab- 
sorbed flux F)y(0.3 — lOkeV) = 5.2x 10“ 14 erg cm -2 s _1 , 
translating into a luminosity of Lx{ 0-3 — 10 keV) = 

1.0 x 10 32 erg s _1 (for a distance of 4.1 kpc), a value 
that is in the range of other young pulsars (Kargaltsev 
& Pavlov 2008). 

The interpretation of the best-fit parameters is dis- 
cussed in Section 4. 

3.3. 7 -Ray Emission 

A 7-ray source is coincident with the posi- 
tion of the PWN with centroid (02000 j ^2000 = 

15 ft, 52 m 45 s , — 56°11 , 46.7"), with 95% confidence radius 
= 3'. 6. The X-ray point source is located near the edge, 
while the PWN is enclosed entirely within the 95% con- 
fidence circle. The origin of the 7-ray emission cannot 
be constrained based on the location of centroid alone, 
since both the pulsar, the PWN, and the SNR itself are 
likely candidates. We discuss these possible scenarios for 
the origin of the 7-ray emission in Section 7. 

The significance of the detection, obtained from the 
evaluation of the peak of the TS map in Figure 5, is 
~ 14.2 <t. The Fermi LAT spectrum from MSH 15-56 is 
shown in Figure 6, along with fits to a power law model 
with an exponential cutoff. The solid curve corresponds 
to the model with the flux in the highest energy fixed 
at the value of the derived upper limit shown in the 
plot, while the dashed curve corresponds to the same 
model with the flux in this bin fixed at a factor of 10 
below the derived upper limit. Assuming a distance of 

4.1 kpc, the range in fit parameters from these two fits is 
E cut « 8.2-43 GeV, T = 1.4- 1.9, andF(> 100 MeV) = 
(1.5 — 3.3) x 10 -8 photons cm' 2 s -1 . The associated lu- 
minosity is Lj(> 100 MeV) = (2.6 — 5.5) x 10 35 erg s _1 . 
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Fig. 5. — Fermi - LAT 7 -ray emission from MSH 15-56 with the 
MOST radio contours overlaid in white. The Fermi source appears 
to be centered at the position of the disrupted PWN. 



Fig. 6 . — Fermi - LAT spectrum of emission from MSH 1556. Also 
shown are model fits for a power law with an exponential cutoff. 
The solid curve is derived by fixing the flux in the highest energy 
bin to be equal to the derived upper limit. For the dashed curve, 
the flux in the bin was fixed at a value ten times below the upper 
limit. 

The ratio L-y/E in this case is somewhat higher than for 
other PWNe, which may in part be explained by the in- 
crease in 7-ray luminosity following the collision between 
the PWN and the reverse shock (Gelfand et al. 2009). 


4. ORIGIN OF THE THERMAL X-RAY EMISSION 

The thermal emission from the northernmost region in 
the XMM image, region labeled (e) in Figure 2a, clearly 
shows silicon and sulfur abundances that are significantly 
higher than solar abundances. This is compelling evi- 
dence that the emission from this region originates at 
least in part from SN ejecta. The best-fit temperature 
and ionization timescale in this region are 0.661 q q 7 keV 
and 1.6 x 10 11 scm” 3 , respectively. The abundances are 
also enhanced in the arc structure in Figure lc, with 
similar values for the best-fit temperature and ionization 
timescale. While the abundances in the other regions 
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are not well constrained, the temperatures and ionization 
timescales are within uncertainties for all regions besides 
region (d) of Figure 2a. The smaller patch of soft X-ray 
emission in region (d) of Figure 2a has a best-fit temper- 
ature of 0.30^0;^ keV, significantly lower than for other 
regions, and a lower limit on the ionization timescale of 
10 13 s cm -3 . 

Based on these spectral properties, it is plausible that 
much of the X-ray thermal emission that we observe with 
Chandra and XMM (all regions besides region (d) in 
Figure 2a) is dominated by ejecta that have been heated 
by the passage of the SN reverse shock, or swept-up or 
mixed-in with the PWN. The enhanced X-ray emission 
at the edges of the PWN likely arises from SN ejecta 
that have been swept up and heated by the nebula. The 
emission in region (d) of Figure 2a may originate from the 
swept-up ISM in the SNR shell. However, there are some 
discrepancies associated with this interpretation that will 
be discussed further in Section 6. 

5. DERIVED SNR PROPERTIES 

We used the best-fit temperature and normalization of 
the thermal component in region (d) of Figure 2a, and 
a distance of 4.1 kpc to calculate the SNR radius R and 
the shock velocity v s , under the assumption of electron- 
ion temperature equilibration, which appears appropri- 
ate for the inferred electron temperature of the shocked 
ISM component (Ghavamian et al. 2007). 


R(pc ) =21 H4.1 

(1) 

v s (km s^ 1 ) = 500 Tg fg keV 

(2) 

Assuming the Sedov (1959) model for the evolution 
of the SNR in the Sedov-Taylor stage, we calculate the 
initial ISM density n 0 , SNR age t, explosion energy E, 
and the swept-up ISM mass M. We find 

no (cm 3 ) =0.1 D,\ /2 

(3) 

t(kyr ) = 16.5 D iA T~^ 2 k eV 

(4) 

E{ 10 51 erg ) = 1.0 D 5 J 2 

(5) 

M(M e ) = 100 DI a 

(6) 


The initial ISM density no was calculated assum- 
ing that ?io = n e /4, where n e is derived from 

the normalization of the thermal spectrum equal to 
10- 14 n e nHV/4ird 2 cm -5 , and a spherical emitting vol- 
ume with a radius of 217 for region (d) in Figure 2a. The 
derived parameters are not unreasonable for a remnant in 
the Sedov-Taylor stage of its evolution. We note that the 
model assumes a uniform ISM density, and a single av- 
erage thermal temperature for the entire swept-up shell. 
The uncertainties in the SNR distance and non-uniform 
ambient density introduce unquantifiable uncertainties 
in the derived parameter values. Based on the derived 
SNR age, we also estimate the transverse velocity of the 
presumed pulsar. The displacement of the point source 
from the geometric center of the shell is approximately 
518, which translates into a transverse velocity of 410 
kms” 1 , within the range of typical pulsar velocities. 


TABLE 3 

DERIVED SNR PROPERTIES 


Property 

Value 

D(kpc) 

4.1 

R(pc) 

21 

ti s (km s” 1 ) 

500 

t(yr) 

16,500 

no (cm 3 ) 

0.1 

£51 ( 10 51 erg) 

1.0 

M( M 0 ) 

100 

V(kn> s *) 

410 


6. PRESENCE OF SN EJECTA 

The enhanced abundances and similar ionization 
timescales of all regions besides region (d) in Figure 2a 
suggest that some of the observed X-ray thermal emission 
originates from SN ejecta that have been either heated 
by the SN reverse shock, or the PWN. The observed ther- 
mal emission is most likely a mixture of emission from 
the swept-up material and ejecta. We can place an upper 
limit on the ejecta mass by assuming that all of this ther- 
mal emission arises from the ejecta (assuming all of the 
ejecta have been shocked). Based on the derived spec- 
tral properties of the thermal emission, we used two dif- 
ferent methods to estimate the mass of this component. 
The two estimates come from two different derivations of 
the electron density n e ; one from the best-fit ionization 
timescale and the derived SNR age, and the other from 
the best-fit normalization of the thermal spectrum. 

Using the age of 16,500 yr and the ionization timescale 
of 2 x 10 n scm” 3 for the presumed SN ejecta component, 
we find the electron density to be n e = 0.38 cm -3 . For a 
spherical emitting volume with a radius of 610 (7 pc for a 
distance of 4.1 kpc) that approximately encompasses the 
X-ray emission visible in the XMM image, we find a total 
hydrogen mass of Mh = 11/ Mg, where / is the volume 
filling factor. If we assume that the ejecta are dominated 
by oxygen, the ejecta density n e j « n e / 8, and the mass 
of each atom is assumed to be 16m p , where m p is the 
mass of a proton. This leads to a total SN ejecta mass 
of M ej = 25 / M 0 . 

We can also derive the electron density n e us- 
ing the best-fit normalization of the thermal spec- 
trum for the same sized region. The fitted nor- 
malization is equal 0.01cm” 5 , which we set equal to 
10~ 14 n e riHV/4:TTd 2 cm” 5 , and find n e / _1 ! 2 = 0.24 cm” 3 
and Mh = 7 f 1 / 2 M 0 . We note that the best-fit pa- 
rameters for this global spectrum are 0.55 keV and 
2 x 10 11 scm” 3 , for the thermal temperature and ioniza- 
tion timescale, respectively, and T = 3.1 for the photon 
index of the non-thermal component that accounts for 
~33% of total unabsorbed flux. 

In order to estimate the total mass of ejecta domi- 
nated by oxygen rather than hydrogen, and correctly 
account for the dependence of the normalization of the 
Bremsstrahlung continuum on the square of the atomic 
number ( Z 2 ), we use the following method. We fit 
the global X-ray spectrum to an absorbed power-law 
model, fixed to the best-fit parameters listed above, plus 
a bremsstrahlung model with a temperature fixed to 0.55 
keV. To reduce the contribution from thermal line emis- 
sion, we limited the spectral fit to energies above 2 keV. 
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Fig. 7. — Top: Model for emission from the evolved PWN in 
MSH 15—50 (see Section 7) along with observed 7 -ray emission 
(and upper limit), radio emission and X-ray flux and spectral index 
from the PWN. Bottom: Time evolution of the PWN magnetic field 
(upper panel) and SNR and PWN radii (lower panel) for MSH 15- 
50. The observed SNR radius is indicated by the horizontal dashed 
line, and the age at which this radius is reached in the model is 
indicated by the vertical dashed line. 

We find a best-fit normalization of ~ 0.01 cm -5 which, 
assuming emission from a metal-rich plasma dominated 
by oxygen, and further assuming n e ~ 8 n,, yields an 
ejecta density n e j ~ 0.02/ -1 / 2 cm -3 and a total ejecta 
mass of M e j « 9.6/ 1/,2 M Q . 

The estimated mass for pure metal ejecta is too high for 
existing progenitor models (Woosley et al. 2002), espe- 
cially considering the fact that the emitting volume used 
in the calculation is only a fraction of the entire SNR 
volume. The enhanced abundances are a clear indica- 
tion that ejecta are present, but the contribution of this 
component may have been overestimated by the fit. A 
lower filling factor for the ejecta would reduce the calcu- 
lated mass. Another possible explanation is that an ad- 
ditional thermal component is contributing to the ther- 
mal spectrum in all regions, but is not accounted for by 
our fits due to the poor statistics of the existing spectra. 
This component may originate from the swept-up ISM 
that contaminates the spectrum along the line of sight. 
On the other hand, multiple components may be present 
if there is macroscopic mixing of metal-rich ejecta with 
hydrogen-rich gas. 

Not accounting for the presence of two separate tlier- 



Fig. 8 . — Models for the broadband emission from MSH 15-5 6 as- 
suming the 7 -ray emission originates from the SNR, dominated ei- 
ther by leptonic process (upper panel) or hadronic processes (lower 
panel). The dotted curves correspond to inverse Compton emis- 
sion, short dashed curves represent nonthermal bremsstrahlung, 
and long dashed curves show the emission from pi 0 decay. 

mal components may be artificially increasing the nor- 
malization of the metal-rich ejecta component that in 
turn leads to a higher inferred ejecta mass. In order 
to test this possibility, we attempted to fit the global 
spectrum of the SNR with a three-component model, 
a power-law, plus two thermal components for the SN 
ejecta and a possible soft ISM component (with a temper- 
ature fixed to 0.3 keV). We find that the presence of this 
additional ISM component could raise the temperature 
of the ejecta component to approximately 0.9 keV, and 
lower its normalization by a factor of ~5, which would in 
turn lower the estimated ejecta mass to ~ 4M 0 . Another 
possibility is that the ejecta are present only in some re- 
gions where we see enhanced abundances, but that most 
of the emission in other regions does not originate from 
ejecta. In this case, it would be difficult to explain why 
both the ejecta and swept-up ISM components have the 
same temperature and ionization timescale. Deeper ob- 
servations of the SNR are required to better constrain 
the spectral parameters and ejecta mass, and confirm 
the origin of the X-ray thermal emission. 

7. ORIGIN OF THE y-RAY EMISSION 

We identified Fermi - LAT emission from 2FGL 
J1552. 8-5609 coincident with MSH 15-5 6 (Figure 5). 
The LAT spectrum is shown in Figure 7, along with the 
radio spectrum (Dickel et al. 2000), the derived flux and 
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spectral index range for the nonthermal X-ray spectrum 
from XMM data, and an upper limit to the TeV emis- 
sion based on HESS non-detection in an accumulated 
~ 40 hours of observations (Hoppe 2008). The evo- 
lution of the composite system was modeled following 
Gelfand et al. (2009), using an injected pulsar spectrum 
with Eq = 3 x 10 38 erg s _1 , using a power law spectrum 
with a e = 2.2 accompanied by a Maxwellian component 
with Ep ea k = 32 GeV comprising 98% of the total elec- 
tron energy. We assumed that 90% of the injected wind 
was in the form of electron/positron pairs, and that the 
PWN expanded into an ejecta core of 25 M 0 driving the 
SNR into an ambient density no = 0.3 cm -3 , roughly 
consistent with that suggested by the X-ray data. 

The resulting evolution of the SNR and PWN radii, 
along with the PWN magnetic field, are shown in Figure 
7. The model reaches the observed SNR radius at an 
age of 1.9 4 yr, by which point the SNR reverse shock has 
already begun significantly compressing the PWN. The 
PWN magnetic field is Bpwn = 34/uG. 

The radio and nonthermal X-ray emission (Figure 7) is 
well-described by synchrotron emission from the evolved 
PWN wind, with the radio emission dominated by the 
Maxwellian electron component, and the X-ray data 
showing the effects of burn-off of the highest energy 
particles due to the increased magnetic field during the 
PWN compression phase. The 7 -ray emission is pro- 
duced through Inverse Compton scattering off of a three- 
component field composed of photons from the CMB, 
starlight (T = 3500 K), and local dust (T = 50 K), with 
equal assumed energy densities. The current spin-down 
power of the pulsar in the model is E = 4 x 10 36 erg s~ : 
with a characteristic age r e = 2.1 x 10 4 yr and an initial 
spin period Pq = 98 ms. We note the limitations of the 
1-D model applied here, particularly given the observed 
asymmetric morphology of the PWN. 

Young pulsars have been shown to form a significant 
component of the Galactic population of GeV 7 -ray 
sources (Abdo et al. 2010), and it is distinctly possible 
that the Fermi LAT emission observed from MSH 15-5 6 
originates directly from the pulsar powering the observed 
radio/X-ray PWN. It is important to note that the ob- 
served cutoff for the power law spectrum is considerably 
higher than that observed for the pulsar population, sug- 
gesting that such an association is unlikely, but pulsation 
searches - both in the gamma-ray band and in other band 
- are of particular importance to address this scenario. A 
search of the ATNF pulsar catalog 10 (Manchester et al. 
2005) reveals no pulsars within more than 2.5 degrees of 
the Fermi LAT source whose spin-down power is larger 
than the observed 7 -ray luminosity at the associated pul- 
sar positions. 

An alternative possibility is that the observed Fermi 
LAT emission originates from relativistic particles accel- 
erated by the SNR itself. We find that models in which 
the emission is produced predominantly by either elec- 
trons, through Inverse Compton scattering of the same 
photon field described above for the PWN scenario, or 
by protons, through p — p collisions and subsequent 7r° 
decay, can both reproduce the observed broadband emis- 
sion (Figure 8 ), though not without some difficulties. 

10 http://www.atnf.csiro.au/research/pulsar/psrcat/ 


Both models require magnetic fields that are much higher 
than expected for the compressed ISM, and also low cut- 
off energies for the particle spectra. For the leptonic sce- 
nario, E etCUt ~ 0.5 TeV and B$nr ~ 80/xG, with nearly 
4 x 10 50 erg in relativistic particles assuming a typical 
electron-proton ratio K e _ p = 10~ 2 . The low maximum 
energy for the electrons is consistent with the high mag- 
netic field, but the field strength itself is quite large. 

For the hadronic scenario, if we assume 40% of the 
available SNR energy (assumed to be 10 51 erg) goes 
into relativistic particles, we require a mean density 
n = 2 cm -2 and K e _ p < 5 x 10 -3 to reproduce the 7 - 
ray spectrum. The associated magnetic field required to 
produce the observed radio emission is Bsnr = 130/zG. 
The proton spectrum requires an exponential cutoff with 
E P}CU t ~ 0.5 TeV, which is surprisingly low given that 
protons do not suffer significant radiative losses. The 
required density is much larger than that inferred from 
X-ray measurements. Such results have been observed 
for SNRs interacting with molecular clouds (Castro & 
Slane 2010), possibly suggesting that those SNRs evolve 
through a clumpy ISM, with a low interclump density 
mitigating the expansion and producing the X-rays, but 
with the relativistic protons that produce the 7 -rays in- 
teracting with both the high density clumps and the 
lower density interclump gas. Such a scenario could also 
apply to MSH 15-56, although there has been no re- 
ported evidence of any interactions with dense clouds for 
this SNR. 

8. CONCLUSIONS 

XMM-Newton and Chandra X-ray observations of the 
composite SNR MSH 15-56 reveal complex structures 
that provide evidence for mixing of the SN ejecta with 
PWN material following a reverse shock interaction. An 
X-ray point source is located at the southern tip of the 
radio PWN, with a trail of hard X-ray emission imme- 
diately behind the point source, indicative of a newly 
formed X-ray PWN. Enhanced X-ray emission traces the 
edges of the radio PWN and may arise from SN ejecta 
that have been swept up and heated by the nebula. 

The X-ray spectra are well described by a two- 
component model, a non-thermal power-law plus a ther- 
mal, non-equilbrium ionization model model with Nh = 
0.5 x 10 22 cm -3 . The thermal emission shows evidence 
for both SN ejecta and the swept-up ISM shell. Se- 
lect regions clearly show Si and S abundances a factor 
of ~3 higher than solar, and a best-fit temperature and 
ionization timescale of 0 . 66 ^q '°7 keV and (1.6 ± 0.2) x 
10 11 s cm -3 , respectively. While the abundances are not 
constrained in all regions, the similarities in temperature 
and ionization timescale suggest that much of the X-ray 
emission originates from SN ejecta that have either been 
heated by the reverse shock or swept up by the PWN. 
We suggest that there also may be a significant under- 
lying ISM component contributing to the spectrum that 
cannot be constrained by the current data. Deeper X-ray 
observations are required to measure the relative contri- 
bution of SN ejecta and the swept-up ISM. We do find 
one region along the southern edge of the SNR shell that 
has a lower temperature of 0.30l°; 04 keV, and appears to 
be in ionization equilibrium. Assuming the Sedov model, 
we derive an SNR age of 16,500 yr. 
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The origin of the observed 7 -ray emission remains 
somewhat unclear. The most likely scenario appears 
to be that it originates from the evolved PWN. How- 
ever, emission directly from the pulsar powering the ob- 
served nebula, as well as emission from the SNR itself 
if the mean magnetic field is very high - and perhaps if 
the remnant is expanding into a clumpy ISM, may also 


provide plausible solutions. Sensitive pulsation searches 
from the pulsar are of considerable importance to ad- 
dress this issue, as are deeper X-ray and molecular line 
measurements of the SNR in order to better assess the 
properties of the ambient medium. 
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